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Abstract 



^ ■ We unravel the structure of the spectrum of the anomalous dimensions of the quark-erluon 

m ; 

CN ! twist-3 operators which are responsible for the multiparton correlations in hadrons and enter as 

Q^ ' a leading contribution to several physical cross sections. The method of analysis is bases on the 

■^ ■ recent finding of a non-trivial integral of motion for the corresponding Hamiltonian problem in 

Ph! multicolour limit which results into exact integrability of the three-particle system. Quasiclassical 

^' expansion is used for solving the problem. We address the chiral-odd sector as a case of study. 
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1 Introduction. 

The studies of power suppressed (in a hard momentum scale) phenomena in physical cross sections 
give an important insight into underlying physics of strong interaction since corresponding quan- 
tities manifest multiparton correlations in hadrons and can shed some light on the complicated 
dynamics of the hadronic substructure. The most promising is the study of twist-three contribu- 
tions in hadronic reactions which enter as a leading effect in certain asymmetries. The planned 
and running experiments will allow to extract these quantities and to face them with theoretical 
predictions. However, the question arises about the scale at which experimental and theoretical 
predictions are confronted. This requires the knowledge of the evolution equations which govern 
the scale dependence of the corresponding correlation functions — a generalization of the familiar 
one-variable parton densities — as well as their solutions, i.e. determination of their eigenvalues 
and eigenstates. While the former problem is quite straightforward since it can be easily resolved 
using standard methods of QCD perturbation theory, provided a basis of operators mixed under 
renormalization is chosen in an appropriate way, and this problem has been tackled by a number 
of authors (see e.g. reviews ||l|, §). On the other hand the second problem is by far more com- 
plicated since it requires the solution of the Faddeev-type many-body equations for the particles 
with pair-wise interaction. Recently a breakthrough has been made, however, by the authors of 
Ref. [^ where they have found an additional non-trivial integral of motion for the problem at 
hand. Even so, when the problem acquires an equivalence to an exactly solvable model the task is 
proved to be too complicated to be solved exactly analytically. Therefore, it is welcome to develop 
an approximate scheme which could allow for a systematical improvement of approximations in- 
volved. The present paper is devoted to the study along this line for the chiral-odd twist-three 
evolution equations, i.e. the ones related to the hadron structure functions /il(x) and e{x) [Q. 

2 Quasi-partonic operators and their evolution. 

The twist-3 sector is exceptional as compared to even higher twists since in the former case contrary 
to the latter we can reduce the renormalization group analysis to the study of the UV divergencies 
of the so-called quasi-partonic operators ||^ which form a compete basis of functions. The twist 
of these objects equals the number of fields which the composite operators are constructed from. 
This means that due to the fact that in leading order only a pair-wise interaction of partons is of 
relevance, the corresponding kernels in sub-channels are the familiar twist-2 non-forward evolution 
kernels. 

The generic form of the evolution equation for twist-3 correlation functions in the momentum 
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fraction formalism looks like 

l/ — Z{xi,X2,x^)=j \{dx,5 [y^Xi-C] K {{xi}\{Xi}) Z{x\,X2,x^ 

with the constraint Xi+ X2+ X'^ = C, for the momentum fractions of partons imposed on the both 
sides of this equation. The variable C, stands for the +-component of the t-channel momentum: 
the limit C = corresponds to the usual DIS kinematics (and Eq. (|l|) to a generalized DGLAP 
equation), and C = 1 to the exclusive one (Brodsky-Lepage equation for baryon distribution 
amplitude). The kernel K in leading order has a pair-wise structure: 



K {{xi]\{x'i}) = ^K (xi, Xj\x'^, x'^ 6 (xi + Xj - x- - x^) 



(2) 
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with K f Xj, Xj\x',i^ x',j\ being the interaction kernel of two nearby particles. 

For the purposes of the present study we need the quark-anti-quark and (anti-)quark-gluon 
kernels with non-contracted Lorentz, Dirac and colour indices. They can be decomposed into 
independent structures as follows 

^i'j';a'P' \Xl, X2|X;^, XgJ 
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for the quark-gluon one. The ellipsis stand for other colour structures which, however, are irrele 
vant for our present consideration. Here as usual the ±-subscripts mean the projection on two op 
posite tangents to the hght cone, and the transversity tensor r^^.^^ = \ [gipQua + ^t^i - ^^^i 
is constructed from the 2D metric tensors of the transverse plane 5^^^ = g^j,y — rip^rvl — riyn* . 
From diagrams in Fig. |]we easily get^ (cf. [0, H, 0) 
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^Here Q^^^ i^{xi,X2, ■■■,Xn) = /_^ ^a™ Hfe^i ("^2;fe — 1 + iO) ''". A detailed discussion of the properties of 
these functions can be found in Ref. 101. 
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Figure 1: The diagrams contributing to the pair- wise quark- ant i- quark kernel at one- loop order in 
the light-cone gauge, B^ = 0. The blobs on the lines stand for the wave function renormalization 
counterterm. 
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Figure 2: Same as in Fig. |l| but for the quark-gluon kernel. The graph with the crossed fermion 
propagator corresponds to the contact-type contribution arising from the use of the Heisenberg 
equation of motion for the quark field. 
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and from Fig. ^ we have (cf. 
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Cf - ^) ^e?i(xi, xi - x',) + ^(/5o + ^Ca)5{x, - x[). (8) 

With these results at hand the construction of the evolution equation for the quasi-partonic 
operators is almost trivial since it requires a mere evaluation of simple traces of Eqs. (^^0) with 
tensor structure of the composite operator under study. 

The correlators related to the chiral-odd structure function e{x) and hilx) are given by light- 
cone Fourier transformation P] 

Z{xr, X2 = -X, - X3, X3) = / ^^e^-i-i+-3X3^;,| 1 ^^{k,, 0, k^) ± Z{-k,, 0, -k,)} \h) (9) 

of the non-local operators 

Z{ki, K2, K3) = -ilj{K3n)af^^ t''gGl^{K2n)-^{Kin). (10) 

The "+" and "— " signs in Eq. (|^) correspond to 1 and 75 structures and, i.e. to e and hi functions, 
respectively. Thus, we have finally the evolution kernel which governs the scale dependence of the 
correlator Z{xi,X2,Xs) 

^°'' {M\{<}) = -^{''^) {xi,X3\x[,x',) 6{x, + X3 - x[ - 4) (11) 



+'^^K^S) (xi, X2\x[, x'2) 5{xi + X2- x^ - x'2 



+'^^Ky^) (Xs, X2IX3, X2) 5{Xs +X2- X3 - X2) - ^S{xi - X[)S{X2 - X2) 
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where we have added the charge renormalization piece due to presence of the coupling constant 
in the definition of the composite operator ([T0|) . Using explicit expressions for pair-wise kernels 
we can find that this is exactly the result of Ref. obtained by a different method. 

3 Consequences of conformal invariance. 

It is well-known that the tree-level conformal invariance of the theory is enough for diagonality of 
the one-loop anomalous dimensions matrix of the conformal operators 0-[|13[- This means that 
the pair-wise kernels displayed above can be diagonalized in the basis of the Jacobi polynomials, 
P) . Namely, due to the support properties of the kernels we can write 



dxidx 



5(xi + X2- x\ - 4 )Pi^'''^''^ (x_^_^2\ <i{M}x^ (xi, xalx' x') (12) 

\Xi-hX2/ 
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where j is a conserved (in leading order) quantum number related to the conformal spin of the com- 
posite two-particle operators, i.e. the eigenvalue of the Casimir operator of the collinear conformal 
algebra su{l,l). The eigenvalues of the kernels are^ 

"%,. = (c^-^)(2^(j + l) + 2^(j + 3)-4^(l)-3), (13) 

'%)^ = (c^-^)(4^(j + 2)-4^(l)-3), (14) 

- (2ih( 



'%)^ = ^i2^PiJ + l) + 2^PiJ + A)-A^Pil)) 



+ (^--^) (, + l)(ri^)(. + 3) +^^°-^^-)-(^^) 
'%)^ = ^(2V^(j + 2) + 2^(j + 3)-4^(l))-(c^-^)^ + i(/?o-3C^). (16) 

where a{j) = (— l)-*. Using this we can deduce the following representation of the kernels |^ 

' 2f^^ CU,(z/i,Z/2) ' \Xi+X2) " ' \X[ + X'J' 

where w{xi,ui\x2,iy2) = Xi'x2' and Wj(z/i,z/2) = ^2jZ[^+u7+i]mj+^tll2+i) ' *^^^ ue = de + sg - 1, 
and de and se being canonical scale dimension and spin of the constituent (pe = {q,g}- 

In the basis of local operators the above eigenfunctions correspond to the conformal operators 
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o,i = u^d^ypt^'"'^ (9+ /9+) 01, (18) 

where 0£ is an arbitrary local field. These operators form an irreducible representation in the space 
of bilinear composite operators of the collinear conformal group with generators J^ ^J~ ^J^ . 
Here the step-up J^ = iV+, the step-down J'~ = |/C_ and the grade J'^ = \{V + A^-+) 

operators, constructed from the momentum, P+, the angular momentum, M. ^, the dilatation, 

P, and the special conformal, /C_, generators, form the sm(1, 1) algebra: [J^,J^]- = ±i7^, 
[J'^,J'^]- = —1J^. The vacuum state is Ojj, [J7'~,(9jj]_ = 0. The eigenvalues of the Casimir 
operator J^ = J^{^J^ ~ 1) ~ J^J^ define the conformal spin, J12, of the state, i.e. [JT"^, Oji\_ = 

Jl2{Jl2 + l)Cj7 with J12 = 3 + \{Vl + V2)- 

Therefore, the evolution equation can be reformulated into the eigenvalue problem for the 
three-particle system in a basis of local operators 

TiIqCD -^ = ^QCD >^, (19) 



^Note that from TV = 1 SUSY Ward identities it follows that ««7'^^ = i^'-i'^ ^ for j =even, "'^T^j+i = ''^7'^^ for 
j =odd at leading order and, therefore, as can be easily seen from the results displayed below in TV = 1 super 
Yang-Mills theory the Hamiltonian of interaction of n-particles with the same helicity is equivalent to the one of 
the XXXs=-i spin chain. 



with Hamiltonian 

Hqcd = '"HfMs) + 'mU-h^) + '"HU-h^) - ^, (20) 



where we have defined "^^yj.i = 2^'^{J) and '^^j_3/2 = 2'^^(J). The operators J's are formally 
determined as solutions of the equation J = J {J + 1). 
In the large A*",, limit Eq. (pOD simplifies into 



^ = Y {^ (•/l2 - ^) + 'A (^12 + ^) + ^ ('^32 - ^) + ^ ('^32 + ^) - 4^(1) -l). (21) 

Thus, the Hamiltonians ( pO| , |2T| ) explicitly manifest the SU{1, 1) invariance of the system. 



4 ^-space. 

Let us explore in full the consequences of the covariance of the problem under the conformal 
transformations. For this we define a space V = {O'^lk = 0, 1, . . . , oo} spanned by the elements^ 

In the representation [i/^'^, x(^)]- = J^'^xi^^) the generators are 

J+ = (, + l), + ,=|, J- = |, J3_l(,^i)^,|, (23) 

with commutation relations: [J^, J^]_ = ±J^, [J'^,J~]_ = —2J^. For a multi-variable function 

x(^i, 02, ■ ■ ■ , On) the operators are defined as J^'^ = Yl^=i Je ' and the quadratic Casimir operator 
is/ = J3(j3_i)_ J+J-. 

Since the spectrum of eigenvalues of the problem (|T9|) have to be real it means that the 



Hamiltonian has to be selfadjoint w.r.t. an appropriate scalar product. We define it as (see [15 
for definition of group invariant measures) 



where 



(X(^l,^2,...,^n)|x(^l,^2,...,^n))= / rfA< x(^l, ^2, • • • , ^n)x(^l, ^2, • • • , ^n), (24) 



dM^H^-^ii-eAT-' and n=\J{n<i] 

e=i 27ri ^^j_ 



^It is worth to note that we work in the space of local operators rather then with the correlation functions 
of composite operator, O, with elementary fields, 0^, {OY[i4>i), or in the language of Refs. ||lj, |^ the so-called 
hat-transformed basis. 



and 6 = 6*. Then the following adjoint properties are obvious 

{x{6,,62,...,6^)\J^''x{0i,6,,...,6^)) = {J^'^x{0i,62,...,6^)\x{6i,62,...,6r,)), 
{X{6,,6,,...,6r,)\nx{0i,62,...,6n)) = {nx{0i,62,...,6nM0i,62,...,6„)). 

Now we are ready to address the question of construction of irreducible representations in 
^-space. 

5 Two-point basis. 

The highest weight vector in the 6'-space depending on two variables, j{^ x(^i) ^2) = 0, is realized 
by translation invariant polynomials^ x(6'i, 6*2) = 612. The descendants are constructed by acting 
with the step-up operator on the vacuum state: (^12) ^12- When transformed to the basis of the 
local operators (E^) they coincide with Eq. ([T8|) up to an overall normalization 






The states with unit norm w.r.t. the scalar product (p^ are 

T{ui)T{u2)T{j + l)r(2j + u, + U2 + l] 



'Pj{0i,62) =n ^{j\ui,U2)6l2, with n'^{j\iy 1,^2) 



r(j + z/i + i)r(j + z/2 + i)r(j + u^ + u2 + i)' 

(25) 
so that {Vj'{6i, 62)\Vj{6i, 62)) = 6j'j. The two-particle Casimir operator 

^2 o2 ^ ^ n // , -,\ ^ / . ^^ d \ ^ iyi + iy2 /'i'i + 1^2 



is obviously diagonal in the basis 

{vAei, e2)\jl2\vAeu O2)) = Sy, [jL],,, with [jj,],, = (j + ^^) fj + ^^ + 1 



(27) 



6 Three-point basis. 

In order to construct an orthonormal basis of three-particles operators we expand them w.r.t. the 

-2 
eigenf unctions of the two-point Casimir operator, say J 12, as follows 



J-j 



k 



Pj;i(^l,^2|^3) = E^fc^3 ' (^y V,{6,,62). (28) 

fc=0 



^We use everywhere the shorthand notation 



From the condition^ J Vj-j{9i, 02|^3) = we deduce the expansion coefficients 

P -N-'(J,\u V ^■) (-imj-j + 1) r(2j + ^, + ^, + 2) 

P,-N (J,jki,z.2,-3)r(^^^)r(^_^._^^^^r(2, + fc + ., + ., + 2)' (2^) 

where the factor 

.r2fj. , r(z/3)r( J - J + i)r(2j + u, + u2 + 2)r(2 j + 1., + ^y, + ^^3 + 2) 

T{J -J + 1^3 + 1)T{J + J + V1 + V2 + 2)r( J + J + z/i + z/2 + z/3 + 2) 

ensures the normahzation of the state to unity (Pj';j'(6'i, 6'2|6'3)|Pj;j(6'i, 6'2|^3)) = Sjij6j'j. From 
the construction it is obvious that J^2 ^^^ total Casimir operator 

J =J^2 + ^23 + -/l3 + E^r^' (31) 

£=1 ^ 

which is the sum of the two-particle Casimir operators in subchannels minus the single-particle 
ones, are diagonal in the basis Vj-ji^Oi, 62103) 

{rj',A0i,e2\93)\jl,\Vj,,{9,,92\93)) = Sj,jS,,,[jI,]j„ (32) 

{Vj>;r{9i,92\93)\f\Vj.,j{9i,92\93)) = 6j>j6j>, [j\, (33) 

where [J ]jj = [J3]jj ([-^3]^^ — l) and [J3]jj = J + |(z^i + z^2 + ^^3 + 3). Since the main quantum 
number J (total conformal spin) is conserved we do not display the dependence on it in matrix 
elements. 

The matrix elements of the remaining generators can be easily evaluated and the result is 

'-2 / '-2 '-2 "2 \ 

{Vj'-j'{9i, 92\93)\J23\Vj.j{9i, 92\93)) = 5j>j [5jij[J23]jj + ^j'j+i[«^23]i+ij + ^i',i-i['^23]j-ij ) (34) 
where the diagonal part is 

[Jl]n = \ [^\ - [^12].. - E ^^ + ^ jj^ (t-^"]- + ^)) ' (35) 

and the non-diagonal elements are 
[4].+i. = }vY^j^(^' + 1)('^ - ^' + ^3), \JW., = '^^|j~)^^ V - J + ^3 + 1), (36) 

-^2 '^2 

with M{J,2) = A^(J, j>i, 1/2, V3)n{2\v\, z/2). So that [J23]i+i,i = [•^23]j,i+i- 

^In what follows operators without subscripts stand for three-particle ones while two-particle operators in 
subchannels are labeled by corresponding indices. 



7 Integral of motion and its quantization. 



On top of the conformal invariance there exists another hidden symmetry of the system (^l]) . The 
beautiful finding made by the authors of Ref. [Q is the identification of an additional conserved 
charge of the three-particle problemQ described by the Hamiltonian ( PT| ) in the limit of large 
number of colours, 

Qt = [jL, jI;V - \ [j\2 + ^Is} , (37) 

which is the hermitian operator and commutes with the Hamiltonian (^) [7i, Qt\- = 0. Every- 
where we have to put z^i = z/3 = 1, z/2 = 2. 

The existence of this additional charge leads to complete integrability of the system (the 
number of integrals of motion equals the number of degrees of freedom). This allows to reduce 
the complicated eigenfunction problem for the Hamiltonian (^Tj) to the more simple one for the 

Qt: 

Q^^ = q^^. (38) 

We will look for the solution in the form of expansion w.r.t. the three-point basis (pSD, i.e. 



^ = T.'^3'Pj;y (39) 

i=o 

Its main advantage comes from the fact that the matrix elements of the charge possess only three 

non-zero diagonals: 

{Vr-AOi, ^2|^3)|Qt|Pj;,(^i, ^21^3)) = h'j (5,',[Qt],, + <5,',,+i[Qt],+i,, + 5,',,-i[Qr],-i,i) , (40) 
with 

[Qt],, = 2 [Ju]n{jl3]n " ^ {[Juh + [jlsh] , (41) 

?2 -, r 4:2 , 9 1 r 42 



[Qt],-+i,, = [[JuWj+1 + [Juh - 2/ [•^23]j+i,i, (42) 

and [Qr]j+i,- = [Qt],o+i- Therefore, the above equation (^) leads to a three-term recursion 
relation for ^j which when solved gives the expansion coefficients as well as quantized values of 
Qt- The former is of the form 

{[Qt],, - qr) ^, + [Qt],,,_i<^,-i + [Qt]j,^,^j+i = 0. (43) 

with the boundary conditions ^^i = !^j+i = 0. The solution for this recursion relation which 
satisfies the boundary conditions exists provided the following constraint is fulfilled 

([QT],,-gT)'<4[Qd,,,+i[QT],,-, (44) 



^See Ref. ||l^ for an exhaustive treatment of the three-quark problem in the context of solution of the Brodsky- 
Lepage evolution equation for baryon distribution amplitudes. 
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Figure 3: The spectrum of the conserved charge qr- Sample trajectories from two sets: for n = 25 
(counts from above, n = 0, 1, . . .) and m = 7 (counts from below, m = 1,2, . . .) are shown by solid 
curves. Long-dashed lines correspond to the analytical formulae (^,^) with n = 2,6. The two 
exact solutions (|5TD separated from the rest of the spectrum by a finite gap are demonstrated by 
short-dashed lines. 



Using the results found so far we can easily deduce from the Eq. ( ^^ the critical eigenvalues 
of the spectrum of the operator Qt for large conformal spin J of the three particle system. 
Substituting the asymptotics for the matrix element of the operator Qt in the three-point basis 



[Qi 



2JV(1 



J— »oo 



2[Q2 



Ji,i+i 



,J~rOO 



J^T^{1 - rf where r = j/J, in Eq. (H) we get < gr < 
r^). The maximum is achieved for Tmax = l/v^- Thus the recurrence relation possesses 
the rising solution for < j ^ 4=J and decreasing one for -j=J<^j < J . We get finally an 
estimate for the quantized values of the Qt for J ^ oo. Namely, the eigenvalues lie in the strip 

1 



< Qt/J < 



(45) 



Note that numerically the upper bound on the spectrum is attained at very high J's. For instance, 
for J = 10^ we have 1.00199 ■ 10"^ < qt/J^ < 0.50531. 

We are now in a position to estimate non-leading in 1/J^ corrections which will generate the fine 
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structure of the spectrum. Let us consider the upper bound of the spectrum. Since the maximal 
eigenvalues are achieved for jmax = 4^^ we consider a small vicinity of jmax: J — 'Jo W ~^ X\/j)- 
We look for the expansion of the eigenf unctions ^j = <?(A) and the charge qt in series 

<PiX) = Y.<P^e)iX)J~'^\ and qrUn) = J^Y.lT\n)J-\ with g^^ = -, (46) 

respectively. Then the difference equation (|43|) is replaced by the sequence of coupled differential 
equations 

<?>;;) (A) + 4 (6 - 4^ - 2A2) cP(o)(A) = 0, (47) 

^'(;)(A) + 4 (6 - 4^ - 2A2) ^(i)(A) = 8A (4v^ - 6 + A^) ^(o)(A), (48) 

The solution of the first one, which satisfies the boundary conditions, i.e. <P(±oo) = 0, is expressed 
by Hermite polynomials 

^(o)(A) = H^ (v^a) e-^^', (49) 

giving the quantized values of (fy 

4\n) = ^U^f2-^-i^, n = 0,l,.... (50) 

Here for a given J, < n < J but the Eq. ( ^Tj) has been derived in the approximation n <^ J and 
thus we can describe the upper part of the spectrum only. The comparison of this approximation 
with the eigenvalues evaluated numerically is shown in Fig. ^ The agreement is reasonable already 
with the first non- leading correction taken into account. We can go to the region n r^ J provided 
a large enough number of terms is kept in Eqs. (^). 

While for low part of the spectrum one can find two exact solutions: 

g--^-l( J) = __ + ( J + 1)2, g-act-2(^) = _^ + ( J + 5)2, (51) 

separated from the rest of the spectrum by a finite gap of order Aq^iJ) oc 0.7 J^ + 40J for large 
J's. 

8 Eigen-energy of three-particle system. 

The eigenfunctions of the integral of motion Qt simultaneously diagonalize the Hamiltonian ( [2T| ) 
and, therefore, the energy of the system is0 



SUn) = ^\2(j:\^A teUmi'-l 



(52) 



^We have used the permutation symmetry of the quark fields so that P^'F = e*''''?', and J^j '^j'^Jdi^Sy ^2|^i) 
e''^Y.j^j'Pj;ji&i,S2\93) since PisT^jijl^i, 021^3) = e^'^T'jy (03,^21^1) where (^ = 0,7r. 
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Figure 4: Same as in Fig. ^ but for the spectrum of energy eigenstates. The two lowest exact 
trajectories are calculated from Eqs. (pSf). 



with e(j) = tp{j + 1) + ip{j + 4) + 2'yE- The upper limit on the spectrum can easily be read from 
the fact that the function ^j is peaked at jmax and thus S"'''''{J) = Nc/2 |2e (j/V^) - §}. But 
this estimate contains non-leading terms as well. 

The energy can be evaluated consistently in WKB approximation as 



£{J,n) 



£^''\J) + J2£^'\n).r 



(53) 



e=i 



Combining Eqs. (|53| ) and (0) with known eigenfunctions allow to determine the coefficients £^^\ 
However, to evaluate next-to- leading WKB correction £^^^ we have to find <P(i) from Eq. (^) as 
well. The solution of this equation is given by 



/oo 
dX'G{X,X')J{X') 
-oo 



(54) 

Here the source i7(A) = 8A f 4^/2 — 6 + A^j ^(o)(A) and the Green function of the homogeneous 
equation is 

I oo g-y2(A2+A'2) 



G(A,A') 



E 



2V2yi -^0 A2™m!(n - m) 



Hr, 



2V2X]H^N2V2X']. 



(55) 



12 



We have omitted, however, in Eq. (0) the solution of the homogeneous equation which might 
be added due to presence of a zero-mode (for m = n), since in the energy the corresponding 



contribution disappears being an odd function of A. Substituting finally the expansion (46) in 
Eq. (|5^) approximating the sum by the integral in the vicinity of the maximum of ^j we find the 
energy to 0{J~') accuracy S = iV,{eW - 3/4 + J-i(2/e(i)<?(o)<?(i) + J e^^^ <P ^o)') / W "^ (o)')} with 
e(°) = 2 In J - In 2 + 27^, e^^) = 2A, e^^) = 4v^ - A^. This gives us finally 

^(o)(J) =41nJ + 47s-21n2--, S^^\n) = 2^2 (3V2 --- nV (56) 

The restriction to this approximation for the energy works somewhat worse then for the g^-charge 
(see Fig. ^ but can be improved routinely. From these results we derive the following expression 
for the energy via the integral of motion at large conformal spin 

£: = In (gr/2) + const + 0( J~^), (57) 

which should be compared with Sxxx = 21ng + const + Oi^J^"^) for the closed XXX quantum 
spin chain [|I^ . 



The lowest trajectories corresponding to the solutions ( |^ reads 

S''''''''-\J) = Y {^^(J + 3) + 27,, - i - j^] , (5J 



-exact-2/ t\ ^c f „ , / t , o\ , r, 1 , ^ 



^exact-2(j) = -^ 2^(J + 3) + 27s - - + 



2 1"' ' '2 J + 3 



And these are exactly the anomalous dimensions found in Ref. [Iq] (see also [IS 



9 Outlook and conclusions. 

To conclude we have constructed a convenient orthonormal three-particle basis in the space of 
local composite operators. Making use of it we have found that the allowed eigenvalues of the 
"hidden" integral of motion Qt and the energy S of the system lie in the strips for large J's: 
J^ ^ (It ^ \J^ and 21nJ<£^<41nJ with the two lowest levels given by Eqs. ( Pl^Sf ) separated 



by finite gaps from the rest of the spectrum. The first non-leading WKB corrections which describe 
the fine structure of the spectra have been found. An improvement of these approximations 
could systematically be achieved by taking into account further O \J~j (^ > 2) terms in the 
recursion relation (^). However, this procedure becomes non-trivial for higher WKB corrections. 
Another fruitful way for determination the energy of the system would be the identification of the 
corresponding three-site integrable open spin chain model [^ with the charge (|37D defined from 
the auxiliary transfer matrix t(A) = Tr {T{\)K~ {X)T^^{—X)K^ [X]) with boundary operators K^ 
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Figure 5: The spectrum of eigenvalues of Hamiltonian (^0[) compared to the range of eigenvalues 
in the multicolour limit (dashed curves). 

satisfying the Sklyanin equation and auxiliary monodromy matrix T(A) = 11^=1-^^ ('^); given by 
the product of Lax operators -^(A) = Al + a^J^ — a^J^ + o'^J^, which fulfills the ordinary Yang- 
Baxter equation. However, since the Hamiltonian 7i ( pT| ) of the system is essentially quantum 
and does not enter into the expansion of t(A) in rapidity A, one has to construct the fundamental 
transfer matrix as well. Once the functional Bethe ansatz is developed it would allow to construct 



the WKB solution ^\, ^ of the corresponding Baxter equation and as a result to deduce higher 
WKB corrections to the Qt in a more economic way together with the corresponding energy. 
Unfortunately the theory of open quantum spin chains [^ is far from being developed to the level 
of periodic models. 

Of course, for real QCD case (finite Nc) the integrability of the system is violated by l/N^ 
corrections in Eq. (120|). These effects manifest themselves in the phenomenon of generation of a 
mass gap between the n = trajectory and the rest of the spectrum (see Fig. ^. 



We would like to thank V. Braun for the talk on the evolution of baryon distribution amplitudes 
p!6| at the Workshop "Structure functions and hadronic wave functions" held in Bad Honnef 
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We are grateful to G. Korchemsky and A. Manashov for very helpful correspondence and the 
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